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ABSTRACT 



We investigate the electron energy distributions (EEDs) and the acceleration 
processes in the jet of Mrk 421 through fitting the spectral energy distributions 
(SEDs) in different active states in the frame of a one-zone synchrotron self- 
Compton (SSC) model. After assuming two possible EEDs formed in different 
acceleration models: the shock accelerated power-law with exponential cut-off 
(PLC) EED and the stochastic turbulence accelerated log-parabolic (LP) EED, 
we fit the observed SEDs of Mrk 421 in both low and giant flare states by using 
the Markov Chain Monte Carlo (MCMC) method which constrains the model 
parameters in a more efficient way. Our calculating results indicate that (1) the 
PLC and LP models give comparably good fits for the SED in low state, but 
the variations of model parameters from low state to flaring can be reasonably 
explained only in the case of the PLC in low state; and (2) the LP model gives 
better fits compared to the PLC model for the SED in flare state, and the intra- 
day/night variability observed at GeV-TeV bands can be accommodated only in 
the LP model. The giant flare may be attributed to the stochastic turbulence 
re-acceleration of the shock accelerated electrons in low state. Therefore, we 
may conclude that shock acceleration is dominant in low state, while stochastic 
turbulence acceleration is dominant in flare state. Moreover, our result shows 
that the extrapolated TeV spectra from the best-fit SEDs from optical through 
GeV with the two EEDs are different. It should be considered in caution when 
such extrapolated TeV spectra are used to constrain extragalactic background 
light (EBL) models. 

Subject headings: Acceleration of particles — Galaxies: active — BL Lacertae objects: 
individual: Mrk 421 — Galaxies: jets — Radiation mechanisms: non-thermal - 
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Gamma rays: galaxies 
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Introduction 



Blazars are the most extreme class of active galactic nucleies (AGNs). Their spectral 
energy distributions (SEDs) are characterized by two distinct bumps. The first bump, which 
is located at the low energy band from radio through UV or X-rays, is generally explained 
by the synchrotron emission from relativistic electrons in a jet that is closely aligned to 
the line of sight. The second bump, which is located at the high-energy band, could be 
produced by inverse Compton (IC 



leptonic model; e.g., 



Bottcher 



scattering of the relativistic electrons (the so-called 
2007f ). The seed photons for IC process can be synchrotron 



photons (synchrotro n self-Compton, SSC; 



radiation fields (EC; 



?lees 



Dermer fc Schlickeiser 



1967 



1993 



Maraschi et al 



Sikora et al 



19921 ) or external 



1994] ). The had ronic mode l 



is an alternative explanation for the high energy emissions from blaza rs (e.g. 



1993 



Miicke et al. 



2003 



Dimitrakoudis et al. 



2012 



Dermer et al. 



Mannheim 



20121 ). 



In the leptonic model, a power-law EED with an ex ponential high-energy 



Tavecchio et al. 



1998 



cutoff (PLC) 



Finke et al. 



or a broken power-law EED is commonly adopted (e.g 
20081 ). The main justification for this EED approximation is that the non-thermal emissions 
from blazars can be described by a power-law spectrum, and the power-law EED can 
be n aturally gener ated in the framework of the shock acceleration (the Fermi I process) 



(e.g. 



Baring 



19971 ) . In recent observations, however, the X-ray spectra of several blazars 
(like Mrk 421, Mrk 501) show significant c urvatu re, which are typically milder than an 



exponential cut-off (IMassaro et al. 



2004a 



20061 ) . Very recently, it has been found that 



gamma-ray emi ssions of many b 



spectrum (e.g., 



Aharonian et al. 



2009 



are success 


hilly 


Abdo et al. 


2010 



itted with a log-parabo lic (LP) 



Ackermann et al. 



20111 ). The LP 



EED is then proposed to model the observed spectral curvature, and such LP EED can be 
generated in the stochastic turbulence acceleration scenario (the Fermi II process) (e.g., 



Tramacere et al. 



2009, 



2011I ). Numerical simulation indicates that stochastic acceleration 
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process may play an imp ortant role in the formation of the particle spectrum in blazar jet 



( IVirtanen fc Vainio 



20051 ) . When the emission mechanisms are determined, the emitting 



EED can be reconstructed from the observed e mission spectra. We can then investiga te the 



acceleration processes acting in blazar jet (e.g., 



Ushio et al. 



2010 



Garson et al. 



20101). 



Blazars are well known for their rapid and large-amplitude variability at all wavebands, 
most prominently at keV and TeV energies. The relations (including correlation and 
time lag) bet ween variabilities at different wavelengths are crucial for constraining jet 



models (e.g., 



Sokolov et al. 



Bottcher & Dermer 



2004; 



Fossati et al. 



2008 



Katarzyhski fc Walczewska 



20101 : 



20101 ) . For instance, if the hard lag is observed, a acceleration process 



could be considered in the model. Different flare patterns indicate different causes of 



(e.g., 


Kirk et al. 


1998 


Chen et al. 


2011 





ge in the injec 


;ion rate and change in the acceleration process, etc 


Kataoka et al. 


2000: 


Graff et al. 


2008; 


Moraitis & Mastichiadis 


2011 



Mr 



t 421 is the closes t known (redshift z =0.031) and the first very high energy (VHE) 



blazar ( IPunch et al 



19921 ). It is classified as high-peaked BL Lac (HBL) according to its 
synchrotron peak location. Mrk 421 is the main target of multi-wavelength monitoring 
campaigns. There are a large number of publications on the multi-wavelength observations 



of this source. Its SED and re 



(e-g- 


Blazeiowski e 


t al. 


2005; 


Albert et al. 


2007: 


2010; 


Bartoli et al. 


2011 


). From these studies, we 



ation of variabiliti e s at different bands are intensively studi ed 



Acciari et al. 



2009 



2011 



Aleksic et al 



with keV flare, however, such relation is very complex. For example, the TeV flux increases 
more than quadratically with respect to the X-ray one, and there is time lag between 
TeV and keV flares. Although SED of Mrk 421 is commonly well fitted by a one- zone 
SSC model, the complex variability behaviors remind us that the realistic model seems 
more complicated. Lately, the multi-wavelength campaign showed evidence for Mrk 421 



-6 - 



in low/quiescent state from 2009 January 19 to 2009 June 1 (lAbdo et al.ll201ll ). During 



this campaign, no significant flare activ ity was seen and t 



the lowest fluxes recorded by MAGIC. 



Abdo et al 



re measured VHE flux is among 



(1201 ll ) therefore claimed that the 



unprecedented complete, 4.5 months average SED observed during this campaign can be 
considered as an excellent proxy for the low/quiescent state SED of Mrk 421. Several 



months later, Mrk 421 was found to undergo i ts one of the 



and gamma-rays bands on February 17, 2010 (jShukla et al 



brightest outburst at X-rays 



the flux correlation at X-ray and TeV band was obse rved ( Shukla et al 



infra - night variability at GeV - TeV band was found (IGalante et al. 



2012 



Raue et al. 



20121 ). 



20121) . During this flaring, 



2011 



2012). and 



Shukla et al 



In this paper we investigate the EEDs and the acceleration processes in the jet of 
Mrk 421 in low state and giant flare state. To achieve our aim, we assume two electron 
distributions, both well motivated by the current particle acceleration models: the shock 
accelerated PLC EED and the stochastic turbulence accelerated LP EED, to model the 
SEDs in the frame of a one-zone SSC model. To more efficiently constrain the model 
parameters and better distinguish between the models, we employ the Markov Chain Monte 
Carlo (MCMC) method instead of a simple ^-minimization procedure to investigate the 
high- dimensional model parameter space systematically. The emission models and MCMC 
method are briefly described in Section [2j In Section [31 we report our results. Finally in 
Section H] we summarize our discussions and conclusions. 
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2. Emission model and MCMC method 

2.1. Emission model 

The one-zone SSC assumes that non-thermal radiation is produced by both the 
synchrotron radiation and SSC process in a spherical blob filled with uniform magnetic field 
(B), which is moving relativistically at a small angle to our line of sight, and the observed 
radiation is strongly boosted by a relativistic Doppler factor Su- The radius of emitting 
blob is R' = <v '™f' pC , where t v ,mm is the minimum variability timescale. Here, quantities in 
the observer's frame are unprimed, and quantities in the comoving frame are primed. Note 
that the magnetic field B is defined in the comoving frame, despite being unprimed. We 



use the methods given in 



Finke et al 



(120081 ) to calculate synchrotron and SSC fluxes. 



The PLC electron distribution formed in the Fermi I acceleration process is 

N'U) = K' el '- S exp(^) for 7min < J < 7max , (1) 

where K' e is the normalization of the EED, s is the electron energy index, and 7^ is 
the electron cut-off energy. 7 min and 7 max is the electron minimum energy and electron 
maximum energy, respectively. In this model, there are eight free parameters, five of them 
specify the electron energy distribution (K' e , 7 min , j' c , 7 max , s) and the other three ones 
describe the global properties of the emitting region (B, R', 8d)- 

The LP electron distribution generated in the Fermi II acceleration process is 

, V ) 7min — 7 — 7c 

N\i) = K' e I v%/ . .v, (2) 

~ c ) Tc — T — 7max ) 



where r is the curvature term of EED (jMassaro et al 



20061 ). In this model, r is an additional 



one besides the eight parameters mentioned above. 
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2.2. MCMC method 



The MCMC method, based on the Bayesian statistics, is well suitable for high 
dimensional parameter space investigation, which is superior to the grid approach with 
a more efficient sampling of the parameter space of interest. The Metropolis-Hastings 
sampling algorithm is adopte d to determin e the jump probability from one point to the 



next in the parameter space (IMackay 



20031 ). The algorithm ensures that the probability 
density functions (PDF) of model parameters can be asymptotically approached with the 
number density of samples. In the follow ing fitting, we will run single chains using the 



Raftery & Lewis convergence diagnostics (IRaftery &: Lewis 



19921 ). and we assume flat priors 



in the model parameters s paces. A brief in troduction to the basic procedure of the MCMC 



sampling can be found in 



Fan et al. 



Neal (1993 ; 



(20 



Gamerman 



oL 



or more details about the MCMC method, 



19971 ); iMackayl (120031 ) . Since the code we used in 



please refer to 
this pa per (ILiu et al. 

and to iLewis fc Bridle! (120021 ) and references therein for a detailed explanation of the code 



2012) is adapted from COSMOMC, we refer the reader to the website 



about the sampling options, convergence criteria and statistical quantities. 



3. Modeling results 



3.1. Modeling the SED in low state 



In this case, we adopt the optical-UV to X-ray data observe d by £w ffi/UVOT/RXT/BAT 



Abdo et al. 



(120 llh . For the Fermi-L AT 



and the Fermi-LAT gamma-rays data reported in 
data, the last two data points are not included in our modeling, due to their very large 
uncertainties. Instead, we use the first two data points measured by MAGIC since the 



x http : / / cosmologist . info/ cosmomc/ 
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extragalactic background light (EBL) absorption on the flux at such energy band is 
negligible. Because the EBL absorption has an effect on the other fluxes measured by 
MAGIC, we also do not take them into account in our modeling. The absorption effect of 
EBL will be discussed in §3.3. The SMA data at 2.3 x 10 11 Hz are used to constrain ^f' min . 
We fix 7^ = 700 for PLC case and 7^ in = 1200 for LP case and 7^ ax = 10 8 for both case. 

For the PLC electron distribution, the one-dimensional (ID) probability distributions 
and two-dimensional (2D) confidence regions (at la and 2a confidence levels) of the model 
parameters, and the SED are shown in Fig. [TJ The fitting parameters are summarized 
in Table [T] with the reduced xl = 1-01 f° r 27 d.o.f. (Note that a relative systematic 
uncertainty of 5 % was added in qua drature to the statistical error of the optical-UV-X-rays 



data reported in 



AbdoetaL 



(120111 )). It can be seen that a very good constraint is derived 
for the spectral energy index s. The parameters j' c , K' e and <5d are well constrained. The 
constraints on B and t v ,mm are not strong, but they are still restricted in relatively small 
ranges. 

For the 2D confidence regions of the parameters, we only show some combinations with 
relatively large correlations. It can be seen that there are negative correlations between 7^ 
and B, B and 8d, as well as t v .mm and 5^, but there are positive correlations between K' e 
and t v ,min and between K' e and 7^. In the following, we try to physically understand the 
correlations among some parameters. In the standard SSC model, the synchrotron SED 
peak flux in the (5-ap proximation is vF y oc N '( , y' r ,)B 2 '~y'J i 5^ and the synchrotron peak energy 



is E s oc JcBdn (e.g 



Tramacere et al 



20111 ). The negative correlation between 7^ and B 



is predicted for the constant E s which can be estimated by t he observed spectra. B ecause 



Tavecchio et al. 



19981 ). where 



B8t> oc v1jv c in the SSC model in the Thomson regime (e.g., 
u s and v c are the peak frequencies of the synchrotron and IC scattering respectively, the 
negative correlation between B and is expected for the observed u s and u c , which 
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implies that the IC scattering is in the Thomson regime. In the SSC model, the radius R' 
of the emitting blob relates to the gamma-ray flux (oc R'~ 2 ) and then can be estimated 
using the observed gamma-ray flux, meanwhile R' ~ c^D^mim hence a negative correlation 
between t v ,min and 5b can be predicted. On the other hand, with B oc E s / (^Sr,), the 
synchrotron peak flux can be rewritten as {vF u ) s oc (K' e /j' c )6 B E 2 , therefore the positive 
correlation between K' e and Y c is expected by the constraint of the observed synchrotron 
peak flux. As to the positive correlation between K' e and t v ,mm, since IC peak flux {vF v ) c 
can be estimated by using the observed gamma-ray flux and {vF v ) c /{uF v ) s oc K'JS&t 



v,min' 



and then Kj/t Vim \ n oc \(vF„) r J(vF„) R ] \v r ]v*\ where t V Tnin B5 D = constant and B6d oc u^/u c 



( iTavecchio et al. 



19981 ) are used. Therefore, the positive correlation between K' e and t v>1 



can be predicted according to the observations. 

For the LP electron distribution, the parameters distributions and SED are shown in 
Fig. [2l The fitting parameters are listed in Tabled] In this case, we have xt = 0-97 (f° r 26 
d.o.f.). The curvature term r is well constrained. It seems that B can be better constrained 
than that in PLC model. The constraints on the other parameters are similar to that 
derived in the PLC model. The explanations on the correlations in the case of PLC still hold 
in the case of LP except for the correlation between B and 7^,. In LP case, the correlation 
between B and Y c becomes weak, which is caused by the effect of r . This correlation 
between B and r is caused by the well-observed spectrum around the synchrotron peak. 

Comparing the results derived in the two EEDs, it can be found that the global 
properties of the emitting blob (B, 5v, R') derived in the two cases are comparable. 
Obviously, we cannot distinguish between the PLC and LP electron distributions in the low 
state based on the above results (see Figs. [1] and [2] as well as Table [p. 
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Fig. 1. — Modeling the SED in low state with the PLC electron distribution. Left: 1- 
D marginalized probability distribution of the parameters. The solid vertical line is the 
expected value, and the 68% limits are depicted by the dashed vertical lines; Middle: 2-D 
confidence contours of the parameters. The contours are for 1 and la levels; Right: the 
best fit to the SED from optical- UV-X-ray to GeV bands. For the low energy component, 
the symbols represent Swift/UVOT (squares), RXT (circles) and BAT (diamonds) data; For 
high energy component, the symbols represent Fermi-LAT (circles) and MAGIC (diamonds) 



data. For the detailed information of the observed data sets, please refer to 
(2011). 



Abdo et al. 
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Fig. 2. — Same as in Fig. CD but with the LP electron distribution. 



Table 1. Fitting parameters for the SED in low state. 



Model 


7'c 


K 


s 


r 


B 


^v,min 


So 


U'JU' B h 




(10 5 ) 


(10 53 /10 42 ) a 






(0.01 G) 


(10 4 s) 


(10) 




PLC model 


2.02±0.11 


1.48±0.39 


2.19±0.02 




3.81±0.56 


3.30±0.75 


3.12±0.32 


34.61 


68% limit 


(1.90 - 2.20) 


(1.04 - 1.91) 


(2.16 - 2.21) 




(3.22 - 4.45) 


(2.50 - 4.18) 


(2.77 - 3.49) 




LP model 


1.47±0.08 


0.93±0.16 


2.43±0.02 


1.68±0.10 


2.05±0.38 


4.90±0.60 


3.34±0.25 


50.71 


68% limit 


(1.37 - 1.60) 


(0.70 - 1.12) 


(2.41 - 2.45) 


(1.58 - 1.78) 


(1.65 - 2.47) 


(4.26 - 5.57) 


(3.07 - 3.62) 





a 10 53 for PLC and 10 42 for LP. 



b Ratio of energy density in electrons to that in magnetic field. 
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3.2. Modeling the SED in flare state 



For the SEP in flaring state, th e Swift/EXT, RXTE/PCA, Fermi-LAT and HAGAR 



data reported in 



Shukla et al 



(120121 ) are adopted. We also fix 7^ in and 7^ ax as we did 
in Section 13.11 The parameters distributions and SED fitting derived with PLC EED are 
shown in Fig.[3j The fitting parameters are listed in Table [2j The resulting xl — 0.79 for 318 
d.o.f. is below unity after a relative systematic uncertainty of 8% was added. The fittings at 
optical and GeV bands are bad. 7^ and s are well constrained, and i V)m i n can be constrained 
in relatively small range. However, the parameters K' e , 5b and B are poorly constrained. 
We cannot obtain the meaningful distribution ranges of the three parameters, thus only the 
68% upper limit are reported in Table [2J Unfortunately, the variability timescale (~ 1 day) 
require d in this PLC model contradicts the observed intra-day variability at GeV - TeV 



bands ( Shukla et al. 



2012 



RaueetaL 



2OI2J). 



Compared to the results in PLC case, it can be seen from Fig. H]that the fittings with 
LP EED is significantly improved with xl = 0-30 (for 317 d.o.f.). Besides, relatively better 
constraints are obtained for all parameters, especially for s, r, 7^ and K' e . Furthermore, the 
LP model with required t Vim i n ~ 8 hours can accommodate the intra-day/night variability 
observed at GeV - TeV bands. It should be noted that the interpretations of the correlations 
of parameters in Section 13.11 still hold here. 



3.3. From low state to giant flare state 



After determining the EED in flare state, we can discuss the EED in low state throu gh 



inves tigating the variations of model parameters with activities. 



Tramacere et al. 



( 2007 



2009|) studied the spectral energy distributions (SEDs) of Mrk 421 in different states in 



the frame of the SSC model and found that there is a negative correlation between the 
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Table 2. Fitting parameters for the SED in giant flare state. 



Model 


J'c 


K 


s 


r 


B 




So 


U'JU'b 




(105) 


(10 53 /10 42 ) 






(0.01 G) 


(10 4 s) 


(10) 




PLC model 


5.82±0.14 


0.04±0.01 


1.65±0.01 




1.07±0.06 


7.98±0.55 


2.99±0.08 


110.08 


68% limit 


(5.66 - 5.97) 


( <0.04) 


(1.64 - 1.66) 




( <H4) 


(7.36 - 8.57) 


(< 3.07) 




LP model 


3.63±0.21 


0.17±0.03 


2.02±0.04 


3.79±0.26 


2.67±0.32 


2.66±0.32 


3.44±0.20 


61.75 


68% limit 


(3.42 - 3.86) 


(0.14 - 0.20) 


(1.99 - 2.06) 


(3.51 - 4.08) 


(2.30 - 3.03) 


(2.30 - 3.02) 


(3.22 - 3.66) 






Fig. 3. — Modeling the SED in giant flare state with the PLC electron distribution. The 
curves are the same as that in Fig. [TJ For low energy component of SED, the symbols denote 
the data from SPOL CCD Imaging/Spectropolarimeter at Steward Observatory (squares), 



the XRT data (circles) and the PCA data (triangles); At high energy component, t 



denote LAT data (circles) and HAGAR data (triangle). Please see 
more details about the data sets. 



r e sym bols 



Shuklaetal 



(120 12[ ) for 
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curvature parameter of the radiatio n spectrum b ~ r/5 and E s , which is expected from the 



stochastic acceleration mechanism. 



Tramacere et al. 



(120 111 ) pointed out that the observed 



negative correlation can be explained by the variation of the momentum diffusion coefficient 
D P {l) = A>o(7/7o) 9 , where q is the turbulence spectrum index (note that a larger value 
of D p implies higher acceleration rate), or the fact that the corresponding EEDs are far 
from the equilibrium where the acceleration dominates over the radiative cooling. They 
also suggested that the curvature increases as the radiative cooling becomes important 
and EED is approaching to the equilibrium during the evolution of EED. From Tables 1 
and 2, if we assume that EEDs in low and giant states both have the LP shapes, we find 
that i c = 1.47 x 10 5 , B = 2.05 x 10" 2 G, 5 D = 34.4, and r = 1.68 in the low state, and 



3.63 x 10 5 , B = 2.67 x 10~ 2 G, S D = 33.4, and r = 3.79 in the giant state. Therefore, 



the ratios of synchrotron peak energies and the curvature parameters in two states are 
[E a ] giant /[E s ] low = [7c 2 ^D]giant/[7 c 2 ^D]iow « 7.7 and [r] giant / [r] low ~ 2.3, which means that 
E s increases with r. This scenario is not compatible with a purely accelerated-dominated 
transition, but the increase in the value of the curvature hints that during the high state, 
the EED is at the equilibrium or very close, and that the cooling is dominating; over the 



acce 



2008 



eration. However, man y other analyses (e.g. 



Becker et al. 



2006 



Stawarz fc Petrosian 



Tramacere et al. 



20111 ) pointed out that in the high state if the EED is close or at the 
equilibrium, the PLC EED would fit the SED better, while our results (Figs. |3] and Hj) show 
that in high state the LP EED fits the SED better. Therefore, the above assumption of 
EEDs in low and giant states both having the LP shapes is not correct. 

Alternatively, we still assume that the EEDs in low and giant states both have the 
LP shapes and the different states are caused by the variation of the momentum diffusion 
coefficient. From Tables CD and [2} it can be found that 7^ and r increases by a factor of ~ 2, 
while 5d and B almost keeps constant from low state to flare state. Since is estimated by 



using the condition t acc (7) = t C ooi(7), where t acc (7) oc p 2 /D p oc 7 2 q /D p0 (ITramacere et al 
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20 111 ) is the acceleration time and i C ooi(7) oc I/7 (where the KN effect of IC scattering is 



neglected) is the radiative cooling time, then we have 7^ oc Dp^ 3 ~ q \ and 7^ oc D p0 in the 
hard-sphere approximation (q = 2). Hence, 7^ increases with D p0 . On the other hand, 
r is inversely proportional to the momentu m diffusion coefficient, i. e., r oc D p $ (e.g., 



Tramacere et al. 



2011 



Massaro et al. 



201 ll ). Therefore, the increase of D p o cannot result in 



the increases of both 7^ and r. 

Since the SED in the giant state using the LP EED can be fitted better in comparison 
with that using PLC EED and the intra-day /night variability observed at GeV - TeV bands 
can be accommodated in LP case (see §3.2 and Tables 1 and 2), therefore the EED in the 
low state may be the PLC shape and the EED in the giant flare state may be LP shape 
here. 



3.4. The EBL absorption 

In this Section, we investigate the EBL absorption. As an example we take the SED 
in low state. Since the EBL absorption becomes important when E >~ 2 TeV for Mrk 421 
(z = 0.031), we compare the TeV spectra predicted by both PLC and LP best-fit models, 
which is shown in the left panle of Fig. [5j It can be seen that the TeV fluxes calculated by 
PLC and LP best-fit models are different when E > 2 TeV. Note that the predicted TeV 
flux is intrinsic flux, therefore, the optical depth for EBL absorption on TeV photons with 
energy E is given by 

r 77 (£) = \n(f int (E)/f ohs (E)) , (3) 

where f mt is the TeV flux calculated by our best-fit model from optical through GeV and 
/ b s is the measured TeV flux. This optical depth can then be compared t o the optical 



depth calculated for the various EBL models (e.g.. iMankuzhiyil et al.ll2010l ). 
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We use the two extrapolated intrinsic TeV fluxes (see Fig. [5]) to calculate t^j. The 
results are shown in right panel of Fig. Obviously, the values of r 77 calculated by using 
the TeV fluxes derived in the LP model is almost two times of those in PLC model when E 
> 2 TeV. In the right panel of Fig [5j for com parison, we also sh ow r 77 calculated by three 



kinds of EBL models: the high level one (e.g . 

— irn i — — 

Franceschini et al.1 12008c iGilmore et al. 



FinkeetaL 



20101 ). the mid dle level one (e.g ., 



Kneiske. fc Dole 



20121 ) and the low level one (e.g. 
20101 ). It can be seen that there are discrepancies not only for the values of r 77 among 
various EBL models but also for those among each EBL model and our results. Therefore, 
when such a method which a intrinsic TeV spectrum is obtained from the extrapolation 
of the best-fit spectrum from optical through GeV is used to constrain the EBL models, 
the emission model of the source should be well determined firstly, at least alternative 
model should be taken into account to compare. Here, due to the large uncertainties of the 
measured TeV data, we cannot obtain meaningful constraints on EBL models or emission 
models. 



As pointed out in 



Massaro et al. 



(120061 ) and 



Tramacere et al. 



( 120091 ). our results 



distinctly show that the curvature in the TeV spectrum can indeed affect the constraints 
on the EBL models. In our fitting, we do not consider the data above ITeV, however as 



discussed in 



Tramacere et al. 



( 120111 ) the property of the curvature in the TeV spectrum is 



more complex. Especially in the extreme KN regime, the curvature of the IC emission is 



large r compared to the Thomson, and it is almost equal to that of the EED ( ITramacere et al. 
201ll ). Therefore, the value of r may be slightly underestimated in the LP case although 
r is mainly constrained by the synchrotron emissions here, and this TeV flux may be 
overestimated. Consequently, a subtle bias on r 77 may be introduced. Therefore, it is 
suggested that in order to constrain the EBL more precisely, the complete TeV data should 
be taken into account in the fitting. In this paper, due to our purpose and the large errors 
of TeV data, our conclusion on the EBL constraint is still robust. 
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Discussion and Conclusions 



Using the MCMC method, we study the SED of Mrk 421 in different states for two 
kinds of EEDs with clear physical meanings: PLC and LP EEDs. Our results indicate that 
the EED in giant flare state is the LP shape and the stochastic turbulence acceleration is 
dominant, while in low activity the EED may be the PLC shape and the shock acceleration 
may play a more important role. This giant flare may be attributed to the re-acceleration 
of these electrons with PLC shape by stochastic process. Basically, we can understand 



Virtanen fc Vainio 



this p rocess from low state to flare state in the scenario proposed by 
(120051 ) that electrons injected at the shock front, then are accelerated at shock by Fermi 
I process and subsequently by the stochastic process in the downstream region. Here, we 
specify this scenario according to our results that the Fermi I process accelerated electrons 
are continuously injected into the emitting blob, and the SSC radiation from the steady 
EED (PLC shape) in the emitting blob is responsible for the emission of Mrk 421 in 
low state; Subsequently, the electrons with PLC shape are re-accelerated by stochastic 
process, and the EED with a significant curvature at high energies (LP) is formed, whose 
radiation contributes to the emission in flare state. A requirement in the scenario is that 
the magnetic field turbulence spectrum must be the case of the so-called "hard-sphere" 
approximation with the spectrum index q = 2 since for cases of Kolmogorov turbulence 
(q = 5/3) and Kraichna (q = 3/2) turbulence, acceleration efficiency depend on the electron 
energy (7'), so that the acceleration efficiency for electrons with 7' ~ 10 5 is very low and 



the escape would be dominant over the acceleration of electrons (e.g. 



Becker et al 



20061 ). 



This scenari o including acce 



flaring (e.g., 



Kirk et al. 



eration process also can account for the spectra hardening in 



1998). This scenario can be examined and the details of the flare 



can be studied in the tim e-dependent model. We will study them in the model includin g 



acceleration process (e.g., 
the coming work. 



Kusunose et al. 



2000 



Katarzyhski et al. 



2006 



Yan et al. 



20121 ) in 
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We notice that the jet of Mrk 421 appears to be p article dominated ( s ee U' IU' U in 



Mankuzhiyil et al 



(120111 ). Finally, 



Table [TP]) . which is consistent with the result derived by 
we stress the caveat on EBL constraints we derived. The best-fit spectrum from optical 
through GeV is often extrapolated into TeV regime, and is considered as the intrinsic TeV 
spectrum. However, our results show that with different emission models, the best-fit SEDs 
from optical through GeV give different TeV spectra. Therefore, it should be considered 
in caution when such extrapolated TeV spectra is used to constrain the EBL models or 
redshift of sourc e. The accurately measured VHE spectrum at E > 2 TeV bands (e.g., 



Raueetal 



20121 ) could put more constraints on EBL models and emission models. 
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Fig. 4. — Modeling the SED in giant flare state with the LP electron distribution. The 
curves are the same as that in Fig. [IJ and the symbols are the same as that in Fig. [3j 
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Fig. 5. — Left: Mrk 421 SED in low state in TeV window; The curves represent the TeV 
spectra calculated by LP model (solid line) and PLC model (dashed line). Right: Opti- 
cal depth r 77 ; The symbols represent r 77 derived by using TeV spectra calculated by LP 
model (di amonds') an d PLC model (cir cles'); The curves r e presen t r 77 calc ulated by EBL 



mode 



of 



Finke et al. 



(120081 ) (dashed), 



Franceschini et al. 



(J2008|) (dotted), 



Gilmore et al. 



fl2012h ( dash-double-do tted for WMAP5+fixed and short-dashed for fiducial WMAP5) and 



Kneiske. fc Dolel (120 10[ ) (dash-dotted). 
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